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ABSTRACT 

The  one  thing  among  many  that  must  be  prepared  for  every  flight  is  the  making 
of  a  flight  plan.  Pilots  must  use  charts  or  graphs  from  the  flight  manual  to  compute  the 
fuel  flow  that  is  essential  to  a  flight  plan.  Since  this  requires  many  steps  of 
interpolation  to  compute  the  specific  conditions  that  cannot  be  read  directly  from  flight 
manual,  it  is  time  consuming  and  increases  the  probability  of  making  a  mistake.  This 
problem  obstructs  the  execution  of  various  mission  changes  and  continuous  sorties. 

A  computer  program  for  personal  computer  or  hand-held  calculator  is  developed 
to  compute  the  desired  fuel  flow  by  modifing  the  equations  for  an  'IDEAL'  airplane. 
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I.  INTRODUCTION 

A.  BACKGROUND 

A  pilot  must  be  furnished  with  flight  performance  information  for  his  aircraft 
that  will  permit  him  to  operate  the  aircraft  efficiently  at  different  airspeeds,  different 
altitudes,  different  values  of  gross  weight  and  different  external  loadings. 

This  information  is  collated  from  a  combination  of  theoretical  determinations 
and  flight  test  data  and  presented  to  the  pilot  in  a  flight  handbook  in  either  a  tabular 
format  or  in  a  series  of  graphs,  or  in  combinations  of  these  methods. 

B.  FUEL  FLOW  RATE 

Probably  the  most  important  parameter  in  optimal  performance  of  an  aircraft  is 
the  rate  at  which  fuel  is  consumed.  It  is  obvious  that  when  the  fuel  is  exhausted  the 
aircraft  will  no  longer  fly,  so  the  pilot  is  very  interested  in  the  rate  at  which  the  fuel  is 
being  consumed.  But  also  of  importance  is  the  fact  that  as  the  fuel  is  burned  the  weight 
of  the  aircraft  is  decreased  and  the  performance  of  the  aircraft  is  changed.  It  is 
therefore  quite  important  that  the  pilot  be  able  to  forecast  the  rate  of  fuel  consumption 
for  each  portion  of  each  flight. 

Fuel  flow  rate  for  a  turbo-jet  aircraft  is  a  direct  function  of  the  thrust  required 
for  the  specified  flight  condition,  and  in  equilibrium  flight  thrust  required  is  equal  to 
the  total  drag  of  the  aircraft. 

C.  TOTAL  DRAG 

The  total  drag  of  an  aircraft  in  equilibrium  flight  may  be  subdivided  into  induced 
drag  (the  drag  associated  with  the  generation  of  lift)  and  parasite  drag  (the  drag 
associated  with  the  resistance  of  the  fuselage  to  motion  through  the  air).  The  induced 
drag  is  a  function  of  the  altitude  of  the  aircraft  (air  density)  and  the  weight  of  the 
aircraft,  and  the  parasite  drag  is  also  a  function  of  the  altitude  as  well  as  being  a 
function  of  the  profile  drag  of  the  aircraft.  The  aircraft  drag,  which  varies  with  the 
configuration  and  external  loading  (weapons,  fuel  tanks,  et  cetera)  is  usually  expressed 
in  terms  of  a  Drag  Index  (DI),  with  the  zero  Drag  Index  usually  being  a  measure  of 
the  drag  of  the  aircraft  in  the  'clean'  configuration. 
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D.  THRUST 

Although  it  is  the  thrust  of  the  engines  that  provides  the  performance  for  a 
turbo-jet  aircraft,  as  long  as  there  is  sufficient  thrust  available  to  provide  the  required 
performance,  the  pilot  is  more  interested  for  the  flight  planning  in  the  fuel  flow  rate. 

E.  FLIGHT  PLANNING 

To  plan  a  flight,  the  pilot  must  determine  what  will  be  the  loading  of  the  aircraft 
(gross  weight)  and  what  effect  this  loading  will  have  on  the  drag  (Drag  Index).  The 
pilot  then  enters  the  flight  handbook  to  determine  the  fuel  flow  rate  for  the  airspeed 
and  altitude  at  which  he  desires  to  fly. 

If  the  gross  weight,  Drag  Index  and  airspeed  are  exactly  the  same  as  values  listed 
in  the  handbook  tables  or  exactly  the  same  as  the  valued  on  the  graphical  lines,  the 
determination  of  the  fuel  flow  rate  is  simply  a  matter  of  reading  (correctly)  the  values 
from  the  information  present  in  the  handbook. 

If,  however,  any  of  the  input  parameters  are  not  identical  to  the  handbook  data 
points,  the  pilot  must  interpolate  between  the  values  presented.  If  weight,  Drag  Index, 
airspeed  and  altitude  are  all  off-presentation  points,  the  pilot  may  be  forced  to  do 
several  cross  interpolations.  This  may  even  involve  interpolations  between  values  listed 
on  different  pages.  Not  only  is  this  a  time  consuming  and  difficult  task,  but  the 
possibility  for  error  is  expanded  with  each  interpolation. 
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II.  APPROACH  TO  THE  PROBLEM 

A.  STATEMENT  OF  THE  PROBLEM 

A  method  was  sought  that  would  permit  the  pilot  to  rapidly  obtain  the  fuel  flow 
rate  for  various  flight  combinations  of  gross  weight,  Drag  Index  and  airspeed.  And, 
regardless  of  the  ease  of  access  to  the  information  and  completeness  of  the  data 
obtained,  this  method  must  provide  fuel  flow  rate  at  the  same  degree  of  accuracy  as  the 
flight  handbook. 

In  addition,  it  was  planned  that  the  method  of  data  retrieval  should  be  as  simple 
as  possible  with  an  ultimate  goal  of  using  a  hand-held  computer  for  the  solutions. 

B.  BASE  LINE  APPROACH 

The  standard  method  of  using  the  flight  handbook  is  for  the  pilot  to  search  and 
follow  several  lines  and  tables  of  the  handbook  to  determine  the  required  fuel  flow  rate. 
This  requires  that  altitude,  gross  weight  and  Drag  Index  must  be  considered,  but  this 
method  is  too  long,  takes  too  much  effort  and  is  open  to  errors. 

The  procedure  can  be  simplified  by  computerizing  the  flight  manual.  Most 
attempts  at  this  method  have  involved  the  development  of  curve-fitting  routines  that 
will  supply  numerical  modeling  of  the  handbook  data.  Although  the  end  use  of  these 
equations  is  quite  rapid,  the  development  of  these  routines  is  a  long  and  arduous  task. 
Inasmuch  as  each  equation  usually  fits  only  a  small  part  of  the  data,  the  computer 
programs  for  such  a  method  are  usually  quite  lengthy. 

The  approach  undertaken  in  this  project  was  the  use  of  the  equations  for  an 
'IDEAL'  aircraft  to  develop  simplified  equations  that  can  be  used  in  a  short  computer 
program. 

Three  steps  were  undertaken  in  the  development  of  these  equations: 

a)  The  fuel  flow  rate  data  in  the  handbook  is  converted  to  thrust  required. 

b)  The  thrust  required  equations  are  adjusted  so  that  they  represent  the  actual 
aircraft. 

c)  The  corrected  thrust  required  equations  are  then  reconverted  to  flow  rate  as  the 
program  output. 
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C.       DETAILED  METHODS 
1.  Concept  of  K«  and  K2 

Determinations  are  made  for  equilibrium,  unaccelerated  conditions  where  drag 
equals  thrust. 

The  basic  drag  equations  for  an  ideal  aircraft  are  used  [Ref.  l:p.  160] 

CD  =  CD0  +  CDi 

CD.  =  CL2/7t  e  AR 

D  =  (l/2)p  V2  S  CD 

All  of  the  principal  items  of  flight  performance  involve  steady  state  flight 
conditions  and  equilibrium  of  the  airplane.  For  the  airplane  to  remain  in  steady  level 
flight,  equilibrium  must  be  obtained  by  a  lift  equal  to  the  airplane  weight.  And  the 
airplane  drag  equals  the  thrust  required  to  maintain  steady  level  flight. 

The  total  drag  of  the  airplane  is  the  sum  of  the  parasite  and  induced  drag. 
Parasite  drag  is  the  sum  of  pressure  and  friction  drag  which  is  due  to  the  basic 
configuration  and,  as  defined,  is  independent  of  lift.  The  equation  of  parasite  drag  can 
be  expressed  as  a  function  of  squared  velocity. 

D0  =  (1/2)CDq  p  S  V2 

Let  Kj  =  (1/2)CD    p  S 


D0  =  K,  V2 

Induced  drag  is  the  undesirable  but  unavoidable  consequence  of  the 
development  of  lift. 

In  the  airfoil,  the  local  lift  vector  is  aligned  perpendicular  the  the  local  relative 
wind  and  hence  is  inclined  behind  the  vertical  by  the  angle  a  .  Consequently,  there  is  a 
component  of  the  local  lift  vector  in  the  direction  of  Vqq  ;  that  is,  there  is  a  drag 
created  by  the  presence  of  down  wash. 

This  drag  is  defined  as  induced  drag,  denoted  by  Dj  [Ref.  l:pp.  152-155] 

Dj  =  (l/2)p  S  V2  CD. 

CD{  =  CL2  /  (*  e  AR> 

Di  =  (CL)2  /(2k  e  AR)p  S  V2 

13 


=  [(2W/p  S  V2)2  /  (2k  e  AR)]  p  S  V2 

=  (4W2  p  S  V2)  /  (p2  S2  V4  7t  e  AR) 

=  [(2W2)/(pSTCeAR)](l/V2) 

If  the  whole  item  of  right  side  except  the  squared  velocity  is  changed  to  some 
constant,  this  equation  can  be  expressed  function  of  squared  velocity. 

Let    K2  =  (2W2)  /  (p  S  k  e  AR) 

Pi  =  K2(l/V2) 

While  the  parasite  drag  predominates  at  high  speed,  induced  drag 
predominates  at  low  speed.  Figure  2.1  illustrates  the  variation  with  speed  of  the 
induced,  parasite,  and  total  drag  for  a  specific  airplane  configuration  in  steady  level 
flight. 

The  curve  of  drag  or  thrust  required  versus  velocity  shows  the  variation  of 
induced,  parasite,  and  total  drag.  Induced  drag  predominates  at  low  speeds  when  the 
airplane  is  operated  at  maximum  lift-drag  ratio,  (L/D)max  the  total  drag  is  at  a 
minimum  and  the  induced  and  parasite  drags  are  equal. 

The  effect  of  a  change  in  weight  on  the  thrust  required  is  illustrated  by  Figure 
2.2.  The  changes  in  the  drag  curves  with  variations  of  airplane  weight,  configuration, 
and  altitude  furnish  insight  for  the  variation  of  range,  endurance,  climb  performance, 
etc.,  with  these  same  items. 

First,  the  primary  effect  of  a  weight  change  is  a  change  in  the  induced  drag  at 
any  given  speed.  Thus,  the  greatest  changes  in  the  curves  of  thrust  and  power  required 
will  take  place  in  the  range  of  low  speed  flight  where  the  induced  effects  predominate. 
The  changes  in  thrust  and  power  required  in  the  range  of  high  speed  flight  are 
relatively  slight  because  parasite  effects  predominate  at  high  speed.  The  induced  effects 
at  high  speed  are  relatively  small  and  changes  in  these  items  produce  a  small  effect  on 
the  total  thrust  required.   [Ref.  2:p.  99] 

Second,  Figure  2.3  illustrates  the  effect  on  the  curve  of  thrust  required  of  a 
change  in  the  equivalent  parasite  area  f  of  the  configuration.  Since  parasite  drag 
predominates  in  the  region  of  high  flight  speed,  a  change  in  f  will  produce  the  greatest 
change  in  thrust  and  power  required  at  high  speed.  Since  parasite  drag  is  relatively 
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Figure  2. 1     Airplane  Thrust  and  Power  Required. 
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EFFECT  OF  WEIGHT  CHANGE 
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Figure  2.2    Effect  of  Weight  on  Thrust  and  Power  Required. 
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small  in  the  region  of  low  speed  flight,  a  change  in  f  will  produce  relatively  small 
changes  in  thrust  required  at  low  speeds.  The  principal  effect  of  a  change  in  equivalent 
parasite  area  of  the  configuration  is  to  change  the  parasite  drag  at  any  given  airspeed. 
[Ref.  2:p.  101] 

Third,  a  change  in  altitude  can  produce  significant  changes  in  the  curves  of 
thrust  required.  The  effect  of  altitude  on  these  curves  provide  a  great  part  of  the 
explanation  of  the  effect  of  altitude  on  range  and  endurance.   [Ref.  2:p.  101] 

Figure  2.4  illustrates  the  effect  of  a  change  in  altitude  in  the  curves  of  thurst 
required  for  a  specific  airplane  configuration  and  gross  weight.  As  long  as 
compressibility  effects  are  negligible,  the  principal  effect  of  increased  altitude  on  the 
curve  of  thrust  required  is  that  specific  aerodynamic  conditions  occur  at  higher  true 
airspeeds. 

As  being  investigated  above,  the  total  drag,  same  as  thrust  required,  is  the 
function  of  squared  velocity  with  different  constant. 

Total  drag  equation  is 

D  =  K,  V2  +  K2/  V2 

where     Ki  =  Cn    p  S 
1  ^o 

K2  =  (2  W2)  /  (p  S  n  e  AR) 

The  first  step  of  this  project  is  to  determine  the  Kj  and  K2  values  to  obtain 
the  thrust  required.    It  can  be  seen  that  the  Kj  is  the  function  of  altitude  and  drag 
index,  K2  is  the  function  of  gross  weight  and  altitude.  So  if  Kj  and  K2  values  can  be 
computed  with  any  flight  condition,  the  thrust  required  will  be  solved  at  any  airspeed. 
2.  Determining  the  Kj 

To  determine  the  Ki,  use  the  minimum  drag  concept. 

At  the  minimum  drag  point,  the  airplane  can  be  flown  maximum  endurance 
and  the  induced  drag  and  parasite  drags  are  the  same. 

The  first  step  for  determine  the  Kj  is  to  investigate  the  thrust  equation.  The 
thrust  equation  can  be  established  with  fuel  flow  for  each  airspeed  and  some  factor  C. 
THRUST    =    FUEL  FLOW  x  C 

Since  the  fuel  flow  of  the  turbojet  powered  airplane  is  proportional  to  the 
airspeed,  the  thrust  required  is  proportional  to  the  fuel  flow. 
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EFFECT  OF  PARASITE  AREA  CHANGE 
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Figure  2.3     Effect  of  Equivalent  Parasite  Area,f,  Thrust  and  Power  Required. 
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EFFECT  OF  ALTITUDE  CHANGE 
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Figure  2.4    Effect  of  Altitude  on  Thrust  and  Power  Required. 
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The  factor  C  is  the  value  of  the  military  power  devided  by  fuel  flow  at  military 
thrust. 

C  =  MILITARY  THRUST  /  FUEL  FLOW  (  at  military  thrust  ) 

The  military  thrust  and  fuel  flow  information  for  any  condition  can  be 
obtained  from  NATOPS  FLIGHT  MANUAL.  With  this  factor  C,  the  thrust  required 
for  each  speed  can  be  computed  by  fuel  flow,came  from  NATOPS  FLIGHT 
MANUAL,  multiplies  factor  C.  Table  1  and  Figures  2.5  and  2.6  show  the  relationship 
and  the  result  between  thrust  and  velocity. 

TABLE  1 
THRUST  REQUIRED  AND  FUEL  FLOW  VERSUS  VELOCITY 

VELOCITY(KTAS)  FUEL  FLOW(LBS/H)     THRUST  REQUTRED(LBS) 

360  7,725  6,558 

400  9,014  7,652 

440  10,623  9,018 

480  12,462  10,579 

520  14,652  12,438 

560  17,377  14,752 

600  21,393  18,160 

MIL  25,680  21,800 

Example)       military  thrust       ;  21800  lb 

fuel  flow  (at  mil)  ;  25680  lbs/hr 
Factor  C  =  21800/25680 
=  0.84891 


Thus 


Thrust  required  =  0.84891  x  Fuel  Flow 


Maximum  endurance  speed  at  specific  gross  weight,  altitude,  and  drag  index 
can  be  figured  out  with  the  graph  of  maximum  endurance  mach  number  in  NATOPS 
MANUAL  Figure  2.7.  This  true  mach  number  needs  to  be  transfered  to  KTAS  to 
substitute  into  the  drag  equation. 

Next  step  is  to  compute  the  minimum  drag. 

Minimum  drag  can  be  computed  from  the  fuel  flow  at  maximum  endurance 
multiplied  by  the  factor  C.  The  maximum  endurance  fuel  flow  can  be  obtained  from 
the  NATOPS  MANUAL  Figure  2.8. 

And  so  the  minimum  drag  is  the  twice  of  the  parasite  drag. 

Dtotal  -  2Di  =  2Do  -  Tmin 
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MAX  ENDURANCE  MACH  NUMBER 
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Figure  2.7     Maximum  Endurance  Mach  Number. 

23 


MAX  ENDURANCE  FUEL  FLOW 


AIRPLAYS  CONFIGURATION 

IMOIVIOUAl  9**6  IHOf  ItS 


REMARKS 


tNCiMtf  II.      7)  J7^-C«_» 

iCao  iTA»o»*oo«r 


BATE.    I10CT0«I»  <*** 

sat*  tAiis.   FLIGHT  TEST 


0 


■'J«L  £«»Of  i 
»'J«L   OtMlITT 

4.*  UA/CAW 

a 

z 

3 

-                               N  *            ■     ■ 

4. 

•  «r 

o 

■  \ 

1 

-.  » 

,^ 1 

O 

*u 

m 

-! 

mS 

a 

'           .          „_             ...  N 

.. 

o 

j 

-j 

j- 

—     aoM" 

.  ! 

Itf 

3      11—* 

a 


2     »  -<si 


nA: 


H       3 

a 


».  vVv« 


— ~%\\V= 


-XWXX: 


— v  t  v  >  v 


CWT 


T 


a 
*      x 


a 
a. 

•  9 


:»5S 


._-.      ^ 


■XWSV- 

HI      2 

4  O 


Figure  2.8     Maximum  Endurance  Fuel  Flow. 
24 


Thus  the  drag  equation  can  be  expressed. 

Tmm  /  2  -  Do  -  (V2)  CDq  p  S  V2 

where  p  =  density  at  given  altiude 

V  =  velocity  at  minimum  drag 
S  =  wing  area  of  specific  airplane 

Cq    can  be  obtained  from  the  above  equation,  because  all  values  except  Cn    are 
o  ^o 

known,  and  the  drag  equation  is  a  function  of  squared  velocity  with  the  known  Cn  . 

l  o 

The  total  value  of  right  side  except  Vz  is  the  desired  constant  Kj. 

3.  Determining  the  K2 

The  procedure  to  compute  K2  is  just  similiar  to  that  for  computing  Ki.  It  is 
convenient  to  use  the  same  minimum  drag  concept.  At  the  minimum  drag  point,  the 
minimum  thrust  required  is  twice  the  induced  drag  amount. 

Use  the  same  maximum  endurance  velocity  and  minimum  thrust  required. 

The  induced  drag  equation  is 

Di=Tmin/2 

=  (2W2)  /  (7C  e  AR  p  S)(l/V2) 

where         w  =  Gross  weight  of  specific  airplane 
s  =  Wing  area  of  specific  airplane 
p  =  Air  density  at  given  altitude 
v  =  Maximum  endurance  velocity 

Since  the  minimum  thrust  required  and  all  items  of  the  right  side  are  known 
value  except  1/tt  e  AR,  the  value  of  1/rc  e  AR  can  be  solved  easily.  After  obtained  the 
1/tt  e  AR  value,  substitute  this  value  into  the  above  induced  drag  equation,  and 
calculated  the  numerical  value  of  the  right  side  except  the  1/V  .  This  numerical  value  is 
just  the  desired  constant  K2  value. 

But  if  the  raw  Ki  and  K2  values  are  used  to  compute  the  thrust  required, 
there  will  exist  a  deviation  of  more  than  10  percent  at  high  speed.  Therefore  one  needs 
to  modify  the  Kj  and  K2  to  fit  the  curve  to  the  FLIGHT  MANUAL  data. 

The  fitting  method  used  in  this  project  is  the  trial  and  error  method.  The 
Table  2  and  Figure  2.9  illustrates  the  difference  between  the  modified  and  raw  value  of 
Kj,  K2. 
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4.  Relationship  between  Kj,  K->  and  Variation  Factors 

Next  consideration  is  to  make  the  equations  of  Kj  and  K9  functions  of  drag 
index,  altitude  and  gross  weight. 

Since  the  Kj  is  the  function  of  drag  index  and  air  density,  Kj  will  be  changed 
by  different  drag  index  and  different  altitude.  Thus  it  is  needed  to  develop  a 
relationship  of  K  ^  with  the  drag  index  and  altitude. 

To  accomplish  this,  it  is  convenient  to  select  some  reference  altitude  and  gross 
weight,  because  the  relationship  of  Kj  and  drag  index  is  independent  of  altitude  and 
gross  weight.  Each  Kj  value  for  every  drag  index  can  be  calculated  with  the  previous 
procedure.  After  calculation  of  all  of  the  Kj  values,  a  graph  of  Kj  versus  drag  index  is 
plotted  to  find  the  relationship. 

In  this  project,  the  altitude  and  gross  weight  are  selected  to  standard  sea  level 
and  40,000  lbs.  The  final  values  and  relationships  are  shown  in  Table  3  and  Figure 
2.10. 

TABLE  3 
DRAG  INDEX  VS  Kl 

Drag  Index  Kj 

0  0.011534 

20  0.013335 

40  0.014807 

60  0.016780 

80  0.018480 

100  0.020527 

120  0.022361 

Fortunately,  as  can  be  seen  in  the  result  of  the  graph,  the  relationship  between 
Ki  and  the  drag  index  is  linear.  So  the  equation  Kj  as  a  function  of  drag  index  is 

K{  =  (0.0115  +  0.0001  x  Drag  Index) 

In  this  equation,  the  most  important  thing  is  the  slope  of  the  graph,  because 
the  starting  value  of  K  j  was  modified  to  fit  the  published  data. 
The  Ki  equation  with  modified  starting  value  of  K^  is 

K{  =  (0.014  +  0.0001  x  Drag  Index) 

This  equation  was  verified  by  calculating  the  drag  at  each  speed.  The  results 
are  shown  in  Tables  4,  5,  6. 
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TABLE  4 
THRUST  REQUIRED  DEVIATION  AT  DI  20 


ALTITUDE 

0     FT 

GROSS  WEIGHT 

40,000  LBS 

D.I 

20 

MILITARY 

THRUST 

21,800  LB 

THRUST 

FUEL  FLOW  x  0.8489 

FUEL  FLOW 

25,680  LB/H 

AIR  SPEED 
(KTS) 

FUEL  FLOW 
(LBS/H) 

THRUST  (LB) 
HAND  BOOK 

THRUST  (LB) 
CALCULATION 

DEVIATION 

360 

7,725 

6,558 

6,456 

-1.5 

400 

9,014 

7,652 

7,741 

1.2 

440 

10,623 

9,018 

9,198 

2.0 

480 

12,462 

10,579 

10,819 

2.3 

520 

14,652 

12,438 

12,600 

1.3 

560 

17,377 

14,752 

14,536 

-1.5 

600 

21,393 

18,160 

16,625 

-8.5 

MAXIMUM  I 

ENDURANCE 

FUEL  FLOW 

5,600       (LBS/H) 

THRUST  (min) 

4,754       (LB) 

AIR  SPEED 

248         (KTAS) 

VARIAT1 

ON  OF   K 

h 

0.014+0.0001x20  =  0.016 

h 

2.0  x  108 

From  Tables  4,  5,  6,  the  following  can  be  observed. 

First,  the  thrust  required,  i.e.,  calculated  drag,  is  almost  same  as  published 
data  at  low  drag  index  range,  but  the  thrust  required  is  much  less  than  published  data 
at  high  drag  index  range. 


30 


TABLE  5 
THRUST  REQUIRED  DEVIATION  AT  DI  60 


ALTITUDE 

0    FT 

GROSS  WEIGHT 

40,000  LBS 

D.I 

60 

MILITARY 

THRUST 

21,800  LB 

THRUST 

FUEL  FLOW  x  0.8925 

FUEL  FLOW 

24,425LB/H 

AIR  SPEED 
(KTS) 

FUEL  FLOW 
(LBS/H) 

THRUST  (LB) 
HAND  BOOK 

THRUST  (LB) 
CALCULATION 

DEVIATION 

360 

9,286 

8,288 

7,935 

-    4.3 

400 

10,980 

9,800 

9,566 

-  2.4 

440 

13,118 

11,708 

11,407 

-  2.6 

480 

15,539 

13,369 

13,449 

-  3.0 

520 

18,662 

16,656 

15,686 

-  5.8 

560. 

23,639 

21,090 

18,115 

-14.1 

600 

MAXIMUM  ] 

ENDURANCE 

FUEL  FLOW 

6,200        (LBS/H) 

THRUST  (min 

) 

5,534        (LB) 

AIR  SPEED 

240          (KTAS) 

variat: 

[ON  OF   K 

Ki 

0.014+0.0001x60  =  0.02 

K2 

2.0  x  108 

Second,  in  all  drag  index  ranges,  the  thrust  required  is  much  less  than 
published  data  at  high  air  speed  range  (approximately  above  Mach  number  0.8). 

To  compensate  for  the  first  problem,  that  is  low  thrust  required  in  high  drag 
index  range,  it  is  needed  to  increase  the  Kj  value,  but  the  thrust  required  value  is 
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TABLE  6 
THRUST  REQUIRED  DEVIATION  AT  DI  80 


ALTITUDE 

0    FT 

GROSS  WEIGHT 

40,000  LBS 

D.I 

80 

MILITARY 

THRUST 

21,800  LB 

THRUST 

FUEL  FLOW  XQ.9165 

FUEL  FLOW 

23'785  LB/H 

AIR  SPEED 
(KTS) 

FUEL  FLOW 
(LBS/H) 

THRUST  (LB) 
HAND  BOOK 

THRUST  (LB) 
CALCULATION 

DEVIATION 

(ID 

360 

10,104 

9,261 

8,674 

-  6.3 

400 

12,094 

11,085 

10,479 

-  5.5 

440 

14,487 

13,278 

12,512 

-  5.8 

480 

17,419 

15,965 

14,763 

-  7.5 

520 

21,477 

19,658 

17,228 

-12.6 

560. 

600 

MAXIMUM  ENDURANCE 

FUEL  FLOW 

6,650         (LBS/H) 

THRUST  (min) 

6,003         (LB) 

AIR  SPEED 

238          (KTAS) 

Kl 

0.014+0.001x80  =  0.022 
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2.0  x  108 

slightly  larger  than  the  published  data.  The  way  to  compensate  in  both  the  low  and 
high  drag  index  range  is  to  increase  the  slope  of  the  K  j  equation  and  to  decrease  the 
starting  value. 
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The  modified  equation  of  K  j  is 

K{  =  (0.0135  +  0.000132  x  Drag  Index) 

And  the  change  of  Kj  equation  can  be  illustrated  in  Figure  2.11. 

The  second  problem,  that  the  thrust  required  is  much  less  than  the  FLIGHT 
MANUAL  data  in  the  speed  range  higher  than  Mach  number  0.8,  can  be  illustrated  as 
the  airplane  reaches  very  high  flight  speeds,  the  drag  rises  in  a  very  rapid  fashion  due 
to  compressibility.  Since  the  generalized  equation  for  parasite  drag  does  not  account 
for  compressibility  effects,  the  actual  drag  rise  is  typified  by  the  dashed  line  in  Figure 
2.12. 

This  phenomenon  can  be  investigated  from  the  viewpoint  of  drag  coefficient. 
The  lift  and  drag  coefficients  vary  with  the  Mach  number,  however,  the  drag  coefficient 
for  a  fixed  value  of  lift  coefficient  does  not  change  until  the  drag-divergence  Mach 
number  is  reached,  at  which  time  the  drag  coefficient  increases  sharply.  If  wind-tunnel 
tests  of  a  model  of  the  airplane  to  be  analyzed  have  been  conducted,  curves  of  drag 
coefficient  against  Mach  number  for  various  values  of  lift  coefficient  will  be  available 
as  shown  in  Figure  2.13. 

Drag  data  above  the  divergence  Mach  number  cannot  be  calculated  with 
precision,  and  to  achieve  any  accuracy  at  all,  wind-tunnel  tests  or  other  empirical  data 
must  be  relied  on.  It  is  noted  that  even  wind-tunnel  data  may  not  be  reliable  in  the 
transonic  range,  and  therefore  performance  guarantees  in  the  transonic  range,  based 
purely  on  wind-tunnel  data  and/or  calculations,  cannot  be  considered  completely 
dependable. 

In  most  instances,  the  lift-drag  relation  may  be  expressed  in  an  analytic  form 
such  as 

CD  =CD     +  CL2/7t  e  AR 

in  which  case  performance  characteristics  may  be  analytically  determined  if  similiar 
analytical  expression  for  power  available  can  be  written. 
The  equation 

CD  =  CD     +  CL2  /  7i  e  AR 

applies  to  most  aircraft  through  the  Mach  number  range  of  0  to  about  4  and  a  C^ 
range  from  0  to  0.6,  for  both  trimmed  and  untrummed  conditions,  however,  Crj    and  e 
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AIRPLANE  DRAG  VS  VELOCITY 
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Figure  2.12    Typical  Airplane  Drag  Curve. 

become  functions  of  Mach  number  above  Mach  numbers  of  about  M  =  0.8.    For 

supersonic  aircraft  the    value  of  Cn     increases  by  a  factor  of  2  to  3  through  the 

o 
transomc  region  and  thereafter  generally  exhibits  a  slight  decrease.   [Ref.  3:pp.  275-278] 

The  span  efficiency  factor  e  shows  a  steady  decrease  from  M  =  0.8  to  about  M  =  2.0 

and  thereafter  remains  almost  constant  at  about  50  to  60%  of  its  low  speed  value. 

From  M  =  0.8  to  M=  1.2  the  parasite  drag  curve  quite  closely  resembles  a  sinJ  curve, 

and  for  approximation  analysis  may  be  represented  by  the  following  equation 
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DRAG  COEFFICIENT  IN  TRANSONIC 
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Figure  2.13     Variation  of  Cd  in  Transonic  region. 
CD     =  CD       +  (CD       "  CD     )  sin3  [(M-0.8)/0.8] 


oo 


om 


oo 


where     Crv        =  maximum  value  of  parasite  drag  coefficient 
om 
Crx        =  minimum  low  speed  value  of  parasite  drag 
oo 

coefficient. 
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The  variation  of  e  may  normally  be  fitted  by  the  polynomial 

e  =  eQ  -  e2  (  M  -  0.8  )2  +  e3  (  M  -  0.8  )3 
where      eQ  =  low  speed  value  of  e 

e2'  e3  =  constants  which  depend  on  the  airplane. 

With  the  above  theory  and  parasite  drag  equation,  it  is  possible  to  make  the 
drag  equation  between  the  Mach  number  M.8  to  M  1.2. 

D0  =  (l/2)pSCDoV2 

at  above  Mach  number  M.8,  where  Mach  Drag  exist, 

D0  =  (l/2)p  S  [CD       +  (CD       -  CD     )sin3  {(M-0.8)/0.8}  ]  V2 
^oo  uom       ^oo 

Kx  =  (l/2)pSCD 
^o 

D0  =  [Kj  +  Kj  I  CD      (CD       -  CD     )  sin3  {(M-0.8)/0.8)  ]  V2 

oo      uom  oo 

=  Kx  [  1  +  (CD       -  CD     )/CD      sin3  {(M-0.8)/0.8}  ]  V2 
uom       ^oo      ^oo 

Since  Kj  and  K2  are  function  of  air  density,  it  must  be  considered  altitude 

correction  factor  of  the  changing  altitude.   [Ref.  4] 

First  consideration  is  the  lapse  rate  and  vertical  temperature  structure  for  the 

standard  atmosphere. 

The  adopted  value  for  the  lapse  rate(-dT/dH  =  a)  in  the  troposphere  is  the 

following  constant 

a)  In  terms  of  standard  geopotential  meters 

a  =  0.0065°C  m"1 


b)  In  terms  of  standard  geopotential  feet  (a) 

a  =  0.0035661 6°  F  ft"1 

c)  In  terms  of  the  c.g.s.  unit  of  H 


a  =  6.628155  x  10"8  °C  cm"2  sec"2 
Accordingly  in  the  troposphere 

T  =  TQ  -  aH  (b) 

where  a,H  and  TQ  are  in  consistent  unit 
Second    consideration    is    the    relationship    between    pressure    and    vertical 
displacement. 
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The  following  relationships  are  to  be  understood  as  expressed  in  any  system  of 
consistent  units. 

The  air  is  assumed  to  obey  the  perfect  gas  law,  which,  in  the  c.g.s  system  of 
units,  may  be  written  as 

p  =  (1/R)(P/T)  (c) 

where 

p   density  of  air,  gm  cm 

p  pressure,  dynes  cm 

T  absolute  temperature,0 K 

R  gas  constant  for  1  gram  of  dry  air,  ergs  gm     {°K} 

And  the  air  is  assumed  also  to  be  in  gydrostatic  equilibrium  and  to  satisfy  the 
differential  equation 

dp  =  -  g  dz  (d) 

where,  in  c.g.s  units, 

P   pressure,  dynes  cm 

p     density  (specific  mass),  gm  cm 

Z  vertical  distance,  cm 

-y 

g  gravitational  acceleration,  cm  sec" 

The    vertical    displacement    is    herein   expressed    in    units    of  geopotential. 
Geopotential  is  defined  in  differential  form  by  the  equation 

G  dH  =  g  dz  (e) 

where 

Z  altitude  measured  positively  upwards  at  a  point 

g  positive  (absolute  numerical)  value  of  the  acceleration 

due  to  gravity  at  the  point 
H  geopotential  at  the  point 

G  dimensional  constant,  the  amount  of  which  determines 
the  magnitude  of  the  unit  of  H  in  terms  of  length  and 
time.(the  dimension  of  G  are  in  units  of  gZ  per  unit 
ofH.) 

Substituting  equation  (b)  in  equation  (c),  equation  (e)  in  equation  (d),  and 
combining  the  results  gives 

dP/P  =  G/R  (-dH/T0  -  aH) 
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=  G/aR  d(T0  -  aH)/(T0  -  aH) 

Intergrating  the  right-hand  member  between  limits  of  0  and  H  gives 

ln(  P/PQ  )  =  G/aR  ln(  TQ  -   aH/TQ  ) 
Let  n  =  G/aR 


Then 


P/P0  =  (T0  -  aH/TQ  )n  =  (  T/T0  )n  (f) 


From  equation  (f),  using  the  consistent  numerical  values  of  G,a,and  R 
n  =  5.2561  (dimensionless) 

The  Figure  2.14  shows  the  relationship  between  temperature,  pressure,  density 
and  altitude. 

With  several  equation  investigated  in  above,  the  useful  relationship  between 
temperature,  pressure  change  and  vertical  displacement. 

For  this  project,  the  useful  relationships  are  written  as 

T/T0    =  1  -  6.875  x  10"6  H 
P/P0    =  (  1  -  6.875  x  10"6  H  )5-2561 
T  =  518.688  x  (  1  -  6.875  x  1(T6  H  ) 
a2  =  1.4  x  1714.87  xT 

Finally  the  relationship  between  density  of  air  and  vertical  displacement  is 

p/p0  =  (P/RT)/(P0/R0T0) 
=  (P/P0)(T0/T) 

And  it  can  be  seen  by  computing  several  cases  that  the  effect  of  changing 
altitude  becomes  larger  than  normal  density  changing  ratio. 

Let       K9  =(p/p0)L2 

Then,  Kj  is  proportional  to  the  Kg  and  K2  is  inverse  proportional  to  the  Kg  . 

As  was  shown  in  previous  equations,  K2  is  proportional  to  the  squared  gross 
weight,  and  the  effect  of  changing  gross  weight  becomes  larger  than  normal  gross 
weight  changing  ratio,  same  as  altitude  change.  Thus,  if  there  was  a  specific  reference 
value  of  K2  for  specific  gross  weight,  the  relationship  between  arbitrary  gross  weight 
and  reference  gross  weight  can  be  made  easily  as 
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K2  =  K2ref(W/Wref)2-7 

With  given  three  specific  conditions.i.e.  altitude,  gross  weight  and  drag  index, 
the  thrust  required  can  be  computed  for  any  airspeed  condition  with  equations  shown 
above.  Kj  and  K2  values  can  be  computed  with  specific  altitude,  gross  weight  and  drag 
index,  and  the  thrust  required  for  any  air  speed  can  be  computed  with  Ki  and  Kj- 

The  final  goal  of  this  project  is  to  compute  the  fuel  flow  rate  for  some  given 
condition,  i.e.  specific  altitude,  gross  weight,  drag  index  and  air  speed.  The  fuel  flow  is 
thrust  required  devided  by  factor  C.  The  factor  C  is  the  value  that  military  power 
devided  by  fuel  flow  of  that  time,  as  was  shown  in  previous  step. 

C  =  MILITARY  THRUST  /  FUEL  FLOW  (at  military  thrust) 
5.  Factor  C 

To  compute  factor  C,  it  is  necessary  to  find  the  relationship  military  thrust 
change  for  different  altitude  and  fuel  flow  change.  This  is  done  in  two  steps. 

a.   Computing  the  military  thrust 

Altitude  is  one  factor  which  strongly  affects  the  performance  of  the  turbojet 
engine.  An  increase  in  altitude  produces  a  decrease  in  density  and  pressure  and,  if 
below  the  tropopause,  a  decrease  in  temperature.  If  a  typical  nonafterburning  turbojet 
engine  is  operated  at  a  constant  RPM  and  true  airspeed,  the  variation  of  thrust  and 
specific  fuel  consumption  with  altitude  can  be  approximated  from  Figure  2.14.  The 
variation  of  density  in  the  standard  atmosphere  is  shown  by  the  values  of  density  ratio 
at  various  altitudes.  Typical  values  of  the  density  ratio  at  specific  altitudes  are  as 
follows 

Altitude,  ft:  Density  ratio 

Sea  level  1.0000 

5,000  0.8617 

10,000  0.7385 

22,000  0.4976 

35,000  0.3099 

40,000  0.2462 

50,000  0.1532 

If  the  fixed  geometry  engine  is  operated  at  a  constant  V  (TAS)  in  subsonic 
flight  and  constant  N  (RPM)  the  inlet  velocity,  inlet  ram,  and  compressor  pressure 
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ratio  are  essentially  constant  with  altitude.  An  increase  in  altitude  then  causes  the 
engine  air  mass  flow  to  decrease  in  a  manner  very  nearly  identical  to  the  altitude 
density  ratio.  Of  course,  this  decrease  in  mass  flow  will  produce  a  significant  effect  on 
the  output  thrust  of  the  engine.  Actually,  the  variation  of  thrust  with  altitude  is  not 
quite  as  severe  as  the  density  variation  because  favorable  decreases  in  temperature 
occur.  The  decrease  in  inlet  air  temperature  will  provide  a  relatively  greater  combustion 
gas  energy  and  allow  a  greater  jet  velocity.  The  increase  in  jet  velocity  somewhat 
offsets  the  decrease  in  mass  flow.  Of  course,  an  increase  in  altitude  provides  lower 
temperatures  below  the  tropopause.  Above  the  tropopause,  no  further  favorable 
decrease  in  temperature  takes  place  so  a  more  rapid  variation  of  thurst  will  take  place. 
The  approximate  variation  of  thrust  with  altitude  is  represented  by  Figure  2.15  and  can 
be  computed  with  the  equation  as  follows 

Tavail  =  Tssl*/°(I/2) 
where       5  =  Pah  /  Pssl 
ff  =  Palt  /  Pssl 

Some  typical  values  at  specific  altitudes  are  as  follows 

Altitude,  ft:         Ratio  of  (Thrust  at  altitude/thrust  at  SSL) 

Sea  level  1.0000 

5,000 0.8960 

10,000  0.8003 

15,000  0.7114 

20,000  0.6296 

25,000  0.5544 

30,000  0.4855 

Since  the  change  in  density  with  altitude  is  quite  rapid  at  low  altitude,  turbojet 
takeoff  performance  will  be  greatly  affected  at  high  altitude.  Also  note  that  the  thrust 
at  35,000  ft.  is  approximately  39  percent  of  the  sea  level  value. 

The  thrust  added  by  the  afterburner  of  a  turbojet  engine  is  not  affected  so 
greatly  by  altitude  as  the  basic  engine  thrust.  The  use  of  afterburner  may  provide  a 
thrust  increase  of  50  percent  at  low  altitude  or  as  much  as  100  percent  at  high  altitude. 
[Ref.  2:pp.  119-121] 
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When  the  inlet  ram  and  compressor  pressure  ratio  is  fixed,  the  principal 
factor  affecting  the  specific  fuel  consumption  is  the  inlet  air  temperature. When  the  inlet 
air  temperature  is  lowered,  a  given  heat  addition  can  provide  relatively  greater  changes 
in  pressure  or  volume.  As  a  result,  a  given  thrust  output  requires  less  fuel  flow  and  the 
specific  fuel  consumption  is  reduced.  While  the  effect  of  altitude  on  specific  fuel 
consumption  does  not  compare  with  the  effect  on  thrust  output,  the  variation  is  large 
enough  to  strongly  influence  a  typical  variation  of  specific  fuel  consumption  with 
altitude.  Generally,  the  specific  fuel  consumption  decreases  steadily  with  altitude  until 
the  tropopause  is  reached  and  the  specific  fuel  consumption  at  this  point  is 
approximately  80  percent  of  the  sea  level  value. 

Above  the  tropopause  the  temperature  is  constant  and  altitude  slightly 
above  the  tropopause  causes  no  further  decrease  in  specific  fuel  consumption.  Actually, 
altitudes  much  above  the  tropopause  bring  about  a  general  deterioration  of  overall 
engine  efficiency  and  the  specific  fuel  consumption  begins  an  increase  with  altitude. 
The  extreme  altitudes  above  the  tropopause  produce  low  combustion  chamber 
pressures,  low  compressor  Reynolds  Numbers,  low  fuel  flow,  etc,  which  are  not 
conducive  to  high  engine  efficiency. 

The  available  military  thrust  for  this  project  with  various  altitude  is  shown 
on  Table  7  and  Figure  2.15. 

TABLE  7 
MILITARY  THRUST  FOR  DIFFERENT  ALTITUDES 

Altitude  (ft)  Military  thrust  (lbs) 

Sea  level  21,800 

5V000  19,533 

10,000  17,446 

15,000  15,508 

20,000  13,725 

25,000  12,086 

30,000  10,584 

The  appropriate  equation  for  the  curve  is  as  follows 
THRUST  =  21,800  -  0.4  x  ALTITUDE 

b.    Computing  the  fuel  flow  at  military  thrust 

There  are  three  effective  factors  required  to  change  the  fuel  flow.  One  is 
altitude,  another  is  drag  index,  and  the  other  is  gross  weight.    To  get  an  appropriate 
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fuel  flow  equation  as  a  function  of  altitude,  drag  index  and  gross  weight,  find  the 
relationship  between  fuel  flow  and  each  effective  factor,  then  combine  the  relationship 
totally. 

(1)  Relationship  between  fuel  flow  and  altitude. 

According  to  the  NATOPS  FLIGHT  MANUAL  data,  the  fuel  flow  is 
decreasing  gradually  with  increasing  the  altitude  with  other  conditions  constant.  To 
find  a  relationship  between  fuel  flow  and  altitude,  specific  drag  index  and  gross  weight 
are  fixed,  to  select  the  available  data.  Table  8  and  Figure  2.16  show  the  relationship 
when  drag  index  and  gross  weight  are  fixed  to  zero  and  40,000  lbs. 

After  curve-fitting,  the  fuel  flow  equation  stated  as  follows 

Fj  =  26,000  -  0.64  x  Altitude  (at  below    8,000  ft) 

Fj  =  20,840  -  0.37  x  (Altitude  -  8,000  )  (at  above    8,000  ft) 

where 

Fj  =    Fuel  flow  considered  with  different  altitude 

(2)  Relationship  between  fuel  flow  and  drag  index. 

According  to  the  NATOPS  FLIGHT  MANUAL  data,  the  fuel  flow  is 
decreased  as  the  drag  index  increased.  As  was  considered  in  step  (1),  a  specific  altitude 
and  gross  weight  will  be  fixed  to  find  the  relationship  between  fuel  flow  and  drag  index. 

Table  9  and  Figure  2.17  show  the  relationship  when  altitude  and  gross 
weight  are  fixed  to  zero  and  40,000  lbs. 

To  find  the  relationship  between  fuel  flow  and  drag  index  with  different 
altitude,  the  computed  fuel  flow  from  step  (1)  must  be  used.  With  above  two  variable 
coefficients,  the  fuel  flow  equation  is 

F2  =  F»  -  16.65  x  Drag  Index  (at  below  Drag  Index  20) 

F2  =  (Fj  +  333.3)  -  22.1  x  Drag  Index         (at  above  Drag  Index  20) 

where 

Fj  =  Fuel  flow  came  from  step  (1) 
F2  =  Fuel  flow  considered  with  different  altitude  and 
drag  index 

(3)  Relationship  between  fuel  flow  and  gross  weight. 

According  to  the  NATOPS  FLIGHT  MANUAL,  the  fuel  flow 
decreased  a   small  amount  with  increasing   gross  weight.   To   find  the  relationship 
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TABLE  8 
FUEL  FLOW  FOR  DIFFERENT  ALTITUDES 

Altitude  (ft)  Fuel  Flow  (lbs/hr) 

Sea  level  26,000 

4,000  23,520 

8,000  20,840 

12,000  19,000 


TABLE  9 
FUEL  FLOW  FOR  DIFFERENT  DRAG  INDEXS 

Drag  Index  Fuel  Flow 

0  26,000 

20  25,680 

40  25,110 

60  24,425 

80  23,785 

100  23,225 

120  22,745 


TABLE  10 
FUEL  FLOW  FOR  DIFFERENT  GROSS  WEIGHTS 

Grossweight  (lbs)  Fuel  Flow  (lbs/hr) 

35,000  26,000 

40,000  26,000 

45,000  25,995 

50,000  25,990 

55,000  25,985 

60,000  25,980 


between  fuel  flow  and  gross  weight,  a  specific  altitude  and  drag  index  must  be  selected 
as  in  the  previous  steps.  Table  10  and  Figure  2.18  show  the  relationship  when  altitude 
and  drag  index  are  fixed  to  standard  sea  level  and  zero. 

It  can  be  seen  from  Figure  2.18  that  the  relationship  between  fuel  flow 
and  gross  weight  is  linear.  The  final  equation  of  fuel  flow  with  different  three  variables, 
i.e.,  altitude,  drag  index  and  gross  weight,  can  be  written  as  follows 

F3  =  F2  -0.002x(Gross  weight-40,000)  (at  heavier  than  40,0001bs) 

F3  =  F2  (at  lighter  than  40,000  lbs) 

where 

F2  =    Fuel  Flow  came  from  step  (2) 
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F-j  =    Fuel  Flow  considered  with  different  altitude, 
drag  index  and  gross  weight 

Now  the  military  thrust  and  fuel  flow  for  any  altitude,  drag  index  and 
gross  weight  can  be  computed  with  the  above  step  (1)  and  (2).  The  desired  factor  C  can 
be  computed  with  the  military  thrust  and  fuel  flow. 

Factor  C  =  Military  thrust/  Fuel  Flow 

To  compute  the  fuel  flow  for  any  various  specific  condition  is  the  final 
goal  of  this  project.  The  fuel  flow  can  be  computed  with  thrust  required,  i.e.  total  drag, 
devided  by  factor  C.  The  total  drag  and  factor  C  had  been  computed  in  previous  step 
already. 

FUEL  FLOW  =  THRUST  REQUIRED  /  FACTOR  C 

With  the  above  equation,  it  is  possible  to  find  the  required  fuel  flow  for 
different  specific  conditions  without  interpolating  the  flight  manual  data. 
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III.  RESULT 

A.       BASE  LINE  CONDITION 

Given  condition 

Gross  weight  40,000  lbs 

Drag  Index  20 

Altitude  Standard  Sea  Level 

1 .  Converting  velocity  versus  fuel  flow  to  velocity  versus  drag 

Military  thrust  (at  sea  level)  21,800  lb 

Fuel  flow    (at  military  thrust)       25,680  lb/hr 

Thrust  =  Fuel  Flow  x  Factor  C 
C  =  Military  thrust  /  Fuel  Flow 

at  military  thrust 

C  =  21,800/25,680  =  0.84891 

The  thrust  required  for  each  velocity  can  be  computed  by  fuel  flow  times 
factor  C(0.84891).   The  Table  11  shows  the  result. 

TABLE  11 
THRUST  REQUIRED  FOR  EACH  SPEED 

VELOCITY(KTAS)  FUEL  FLOW(LBS/H)     THRUST  REQUIRED(LBS) 

360  7,725  6,558 

400  9,014  7,652 

440  10,623  9,018 

480  12,462  10,579 

520  14,652  12,438 

560  17377  14,752 

600  21,393  18,160 

MIL  25,680  21,800 


2.  Computing  the  minimum  thrust  and  maximum  endurance  speed. 

Maximum  endurance  speed  can  be  obtained  from  maximum  endurance  Mach 
number  graph  in  NAPTOPS  FLIGHT  MANUAL  (see  Figure  2.7).  The  maximum 
endurance  Mach  number  is  M  0.378.    And  maximum  endurance  fuel  flow  can  be 
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obtained  from  the  maximum  endurance  fuel  flow  graph  in  NATOPS  MANUAL  (see 
Figure  2.8). 

The  minimum  thrust  can  be  obtained  from  the  maximum  endurance  fuel  flow 
by  multiplying  with  the  factor  C  computed  in  step  1. 

maximum  endurance  fuel  flow  ;  5,600  lb/hr 
Tmin  =  5'600  x  °-84891  "  4,753.9  lb 

3.  Computing  the  Kj 

At  maximum  endurance  speed,  the  thrust  required  is  equal  to  the  total  drag, 
that  is  twice  the  induced  drag  or  parasite  drag,  i.e.  induced  drag  and  parasite  drag  are 
the  same  at  this  point. 

Tmin  =  Dtotal  =  2Di  =  2Do 
D0  =  CD  (l/2)pV2A 

where       DQ  =  4,753.9  /  2  =  2,376.9  lb 
p  =  0.0023769  lb  sec2/ft4 
V  =  0.378  x  1,116.4  ft/sec 
A  =  530  ft2 

CD0  =  2Do  /  (Pv2a) 
•    °  =  (4,753.9)/(0.0023769  x  422.1 82  x  530) 

=  0.02117183 

D0  =  (0.021 17183)(1/2)(0.0023769)(530)V2 
=  0.01333  V2 
Thus   Kj  =  0.0133333 

4.  Computing  the  K2 

Use  the  same  concept  as  for  Kj 

Tmin  -  Dtotal  =  2Di  =  2Do 
D0  =  2W2  /  (pAV2)(l/7reAR) 

1/rceAR  =  2W2  /  (pAV2Di) 

=  (0.0023769)(530)(422.18)2(2376.9)/(2  x  400002) 
=  0.16677964 
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Dj  =  (2W2/pA)(l/7teAR)(l/V2) 

=  (2  x  40,0002  /  0.0023769  x  530)(0.16677964)(1/V2) 

Thus 

K2  =  4.2365  x  108 
Dtotal  =  Thrust  required  =  K{V2  +  K2/V2 

These  Kj  and  K2  values  are  modified  a  little  to  fit  the  curve  of  hand  book 
data  by  the  trial  and  error  method  introduced  in  Chapter  2.  Thus  the  thrust  required 
equation  can  be  written  as 

D  =  0.01614V2  +  2.0xl08/V2 

But  if  there  were  no  consideration  about  the  Mach  effects,  the  drag  at  high 
Mach  number,  i.e.  greater  than  M  0.8,  will  be  much  less  than  the  hand  book  data.  Cq 
becomes  a  function  of  Mach  number  above  Mach  numbers  of  about  M  0.8.  In  this 
region,  the  parasite  drag  curve  quite  closely  resembles  a  sin  curve,  and  for 
approximation  analysis  may  be  represented  by  the  equation 

cn    =  cn      +  (cn      -  cn     ^ sin3  ttM  -  0.8)/0.8} 
^o  ^oo  uom       ^oo^ 

Thus  the  drag  equation  between  the  Mach  number  M  0.8  and  M  1.2 

D0  =  1/2  pSCDqV2 

D_  =  1/2  pS  [CD       +  (CD       -  CD     )sin3{(M  -  0.8)/0.8}]  V2 
u  ^oo  ^om       ^oo 

K{  =  l/2pSCD 

D0  =  [K,  +  K,/CD      (CD       -  CD     )  sin3{(M  -  0.8)/0.8}]  V2 

u  l  l        oo      uom  oo 

=  K,[1  +  (CD       -  CD     )/CD      sin3{(M  -  0.8)/0.8}]  V2 
1  uom       ^oo      ^oo 

where  Cn       -  Cn       =  0.55 
uom       ^oo 

Cn       =  0.018 
^oo 

D0  =  Kj[l  +  0.55/0.018  sin3  {(M  -  0.8)/0.8}]V2 
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5.  Converting  to  fuel  flow  from  drag 

Fuel  flow  for  each  velocity  can  be  computed  from  thrust  required,  i.e.    total 
drag,  devided  by  factor  C. 

C  =  MILITARY  THRUST  /  FUEL  FLOW  (at  military  thrust) 

Military  thrust  and  fuel  flow  can  be  computed  with  the  equation  that  found  in 
Chapter  2. 

Military  thrust  =  21,800  -  0.4  x  Altitude 
=  21,800-0.4x0 
=  21,800 
Fuel  flow; 

Fj  =  26,000  -  0.64  x  Altitude 
=  26,000  -  0.64  x  0 
=  26,000 
F2  =  Fj  -  16.65  x  Drag  Index 
=  26,000-  16.65x20 
=  25,667 
F3  =  F2  -  0.002  x  (Gross  weight  -  40,000) 
=  26,667  -  0.002  x  (40,000  -  40,000) 
=  25,667 

Thus  C  =  21,800/25,667 
=  0.84934 

Fuel  flow  =  Thrust  required  /  0.84934 

Table  12  and  Figure  3.1  show  the  results  of  the  calculation  and  deviation  from 
hand  book  data. 

B.       VARIATION  NUMBER  1  -  GROSS  WEIGHT 

Given  condition 

Gross  weight       50,000  lbs 
Drag  Index  20 

Altitude         Standard  Sea  Level 
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The  gross  weight  is  changed  only  from  the  base  line  condition.  It  is  necessary  to 
consider  the  change  of  Kj  and  K2. 

K,  =  1/2  pCn  S 

1  ^0 

As  you  can  see  from  the  above  equation,  Kj  is  independent  of  gross  weight. 
Thus  there  was  no  change  in  Ki. 

But  K2  is  a  function  of  gross  weight  (as  discussed  in  chapter  3). 

K2  =  2W2/(p  S  re  e  AR) 

It  can  be  seen  that  the  K2  is  proportional  to  the  squared  gross  weight,  but  the 
effects  of  varing  gross  weight  become  more  largely. 

Thus       K,  =  K9  (50,000/40,000)2-7 

=  2.0xl08x(1.25)2'7 

=  3.653  x  108 

D  =  0.0164V2  +  3.653xl08/V2 

Military  thrust  =  21,8001b 

Fuel  flow; 

Fj  =  26,000  -  0.64  x  Altitude 

=  26,000 
F2  =  Fj  -  16.65  x  Drag  index 
=  26,000  -  16.65  x  20 
=  25,667 
F3  =  F2  -  0.002  x  (Gross  weight  -  40,000) 
=  25,667  -  0.002  x  (50,000  -  40,000) 
=  25,647 
C  =  21,800/25,647 
=  0.85 
Fuel  flow  =  Thrust  required  /  0.85 

Table  13  and  Figure  3.2  show  the  results  of  calculation  and  deviation  from  hand 
book  data. 
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C.      VARIATION  NUMBER  2  -  DRAG  INDEX 

Given  condition 

Gross  weight      40,000  lbs 

Drag  Index  40 

Altitude         Standard  Sea  Level 

If  the  drag  index  was  changed  from  base  line  condition,  what  will  be  effect  on  the 
result?    It  is  necessary  to  consider  the  Kj  and  K2  equation. 

K{  =   l/2PSCDo 

K2  =  2W2  /  (pSrceAR) 

As  can  be  seen,  K2  is  a  function  of  gross  weight,  altitude,  wing  area  and  aspect 
ratio,  not  of  drag  index,  thus  K2  was  unchanged.  But  Kj  is  a  function  of  altitude,  wing 
area  and  drag  index.  Thus  the  appropriate  Kj  value  can  be  computed  with  equation 

K,  =  K,  x(Ki  /K,  ) 

*DI40  ]DI20  ]DI40      *DI20 

Now  it  is  necessary  to  find  the  relationship  of  Kj  between  drag  index  20  and 

drag  index  40.  Furthermore,  if  the  relationship  of  Kj  between  any  drag  index  and 

reference  drag  index,  the  various  Kj  for  different  drag  index  can  be  computed.  To 

figure  out  the  relationship  of  Kj  for  each  drag  index,  let  us  compute  each  Kj  value. 

The  computing  method  is  the  same  as  introduced  in  the  base  line  condition  step. 

T         =  2D 

1min       LtJo 

D0  =  1/2  pSCDqV2 

CD     =  2D0/(pSV2) 

K{  =  1/2pSCDq 

For  example 
when  DI  =  0 

Do  =  Tmin/2  =  4,443.8/2=2,221.9 

CD     =  2(2,22 1.9)/(0.0023769x530x438.92) 
°  =0.018132 


Thus 


Kj  =  l/2pSCD 
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=   1/2(0.0023769x530x0.018132) 
=  0.011534 

Table  14  and  Figure  3.3  show  the  result  and  the  relationship  of  Kj  for  different 
drag  index. 

TABLE  14 
Kl  VALUE  VERSUS  DRAG   INDEX 


DI 

C 

MAXIMUM  ENDURANCE 

THRUST (LB) 
MINIMUM 

CDo 

Kl 

FUEL  FLOW 

KTAS 

0 

0.8385 

5,300 

260 

4,444 

0.018 

0.012 

20 

0.8489 

5,600 

250 

4,754 

0.021 

0.013 

40 

0.8682 

5,920 

247 

5,140 

0.024 

0.015 

60 

0.8925 

6,250 

241 

5,578 

0.027 

0.017 

80 

0.9165 

6,520 

238 

5,976 

0.029 

0.018 

100 

0.9386 

6,800 

233 

6,383 

0.033 

0.021 

120 

0.9584 

7,050 

230 

6,757 

0.036 

0.022 

The  Ki  of  the  base  line  condition  was  modified  to  fit  the  drag  index.  As  you  can 
see  in  Figure  3.3  and  the  investigation  in  Chapter  2,  the  relationship  between  Kj  and 
drag  index  is  linear  and  the  equation  of  Kj,  function  of  drag  index  is 


Ki  =  0.0135  +  0.000132  x  Drag  Index 


Thus 


K 


DI40 


=  0.0135  +  0.000132x40 


=  0.01878 

Thrust  required  =  0.01878V2  +  2xl08/V2 
Military  thrust  =  21,800  lb 


63 


X 

w 

Q 

O 
< 

Q 
> 


o 

"3 


C3 


3 
1/5 


o 


'■J 

3 
CO 


220'0       0200 


9T0"0       9T0"0       tTO'O 

amvA  tx 


SIO'O   OTO'O 


64 


Fuel  flow; 

Fj  =  26,000 

F2  =  (Fj  +  333.3)  -  22.1  x  Drag  Index 
=  (26,000  +  333.3)-  22.1  x  40 
=  25,449.3 

F3  =  F2  -  0.002  x  (40,000  -  40,000) 

=  25,449.3 

C  =  21,800/25,449.3 
=  0.8566 

Fuel  flow  =  Thrust  required/0.8566 

Table  15  and  Figure  3.4  show  the  results  of  calculation  and  deviation  from  hand 
book  data. 

D.       VARIATION  NUMBER  3  -  ALTITUDE 

Given  condition 

Gross  weight     40,000  lbs 
Drag  Index  20 

Altitude  4,000  ft 

If  the  altitude  will  be  changed  only  from  the  base  line  condition,  how  will  it  affect 
the  Kj  and  K2? 

The  equations  of  Kj  and  K2  are 

Kj  =  l/2PSCDo 

K2  =  2W2/(pSrceAR) 

As  can  be  seen  in  the  above  equation,  the  Kj  and  K2  both  are  functions  of  the 
air  density,  and  air  density  is  a  function  of  altitude.  Thus  the  Kj  and  K2  must  be 
changed  with  varying  altitude. 

The  useful  relationships  came  from  the  ICAO  Report  written  as 

p  =  P/RT 

T/T0  =  1  -  6.875xl0-6H 

P/P0  =  (1  -  6.875xl0-6H)5-2561 
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T  =  518.688  x  (1  -  6.875xl0'6H) 
And  the  relationship  between  density  of  air  and  vertical  displacement  is 

p/p0  =  (P/RT)(P0/RT0) 

=  (P/P0)(TQ/T) 

Since  the  effects  of  changing  altitude  become  larger,  it  is  necessary  to  modify  the 
equation. 

Let       K9  =  (p/p0)L2 

Then  Kj  is  proportional  to  the  Kg  and  K2  is  inverse  proportional  to  the  Kg. 
Thus  in  this  case 

p/p0  =  (P/P0)(T0/T) 

=  (  1  -  6.875xl(r6x4,000)5-2561/(l  -  6.875xl0'6  x  4,000) 
=  0.88809 

K,  =  Ki       x(0.88809)L2 
ssl 
=  .0139976 

,1.2 


K2  =  K2      /  (0.88809)J 
=  2.306xl08 
Thrust  required  =  0.0139976V2  +  2.306xl08/V2 
The  military  thrust  will  be  reduced  as  discussed  in  Chapter  2. 

Tavail=Tssl5/°(1/2) 

5  =  P4,000ft  /  Pssl 

=  1,827.69/2,116.22 
=  0.8636578 

a  =  P4,000ft  /  Pssl 

=  0.0021109/0.0023769 
=  0.8880895 

Tavail  =  21,800  xO.8636578/0.88808951/2 
=  19,779 
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Fuel  flow; 

Fj  =  26,000  -  0.64  x  Altitude 

=  26,000  -  0.64  x  4,000 

=  23,440 
F2  =  Fi  -  16.65  x  Drag  Index 

=  23,440  -  16.65  x  20 

=  23,107 
F3  =  F2  -  0.002  x  (Gross  weight  -  40,000) 

=  23,107  -  0.002  x  (40,000  -  40,000) 

=  0.86463 

C  =  19,979/23,107 
=  0.86463 

Fuel  flow  =  Thrust  required/0.86463 

Table  16  and  Figure  3.5  show  the  results  of  calculation  and  deviation  from  hand 
book  data. 

E.       VARIATION  NUMBER  4  -  WEIGHT,  DRAG  INDEX  AND  ALTITUDE 

Given  condition 

Gross  weight      50,000  lbs 
Drag  Index  40 

Altitude  8,000  ft 

If  ?U  conditions  are  changed  from  base  line  conditions,  it  is  neceassary  to 
consider  for  each  changed  condition  step  by  step,  same  as  previous  sections. 

At  first,  think  about  gross  weight  change.  The  changing  gross  weight  affects  the 
^2  only,  because  Kj  is  independent  of  gross  weight. 

K2  -  2W2  /  (pSrceAR) 

As  you  can  see  in  section  3.2 

K-,  =  K-,  (50,000/40,000)2J 

250,0001bs  240,0001bs 

=  2.0xl08(1.25)2J 

=  3.653xl08 
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Next,  think  about  the  drag  index  change.  If  the  drag  index  is  changed,  it  is 
necessary  to  consider  Ki  only,  because  the  K2  is  independent  of  the  drag  index  change. 
K  J  equation,  function  of  drag  index  is 

Kj  =  0.0135  +  0.000132  x  Drag  Index 

Thus 

K,     =  0.0135  +  0.000132x40 
^140 

=  0.01878 

And  finally,  think  about  the  altitude  change. 

As  you  can  see  in  section  3.4,  the  altitude  change  affects  both  Kj  and  K2, 
because  Kj  and  K2  are  function  of  air  density.  The  equation  in  section  3.4  can  be  used 
in  this  step. 

p/p0  =  (P/RT)(P0/RT0) 

=  (P/P0)(T0/T) 


Let   K9  =  (p/p0)12 

K,  =  f(p) 
K2  =  f  (1/p) 
Thus 


Kl  _  KlDI40x(p8,000ft/Pssl) 


1.2 

=  0.01878  x  (0.0018683/0.0023769)1-2 
=  0.01406755 

K2  =  K2cnnnnlK   /(Ps.OOOft/Pssl)1' 


-50,0001bs 
653xl08  ( 

=  2.9157xl08 


=  3.653xl08  /  (0.0018683/0.0023769)1'2 


Thrust  required  =  0.01407V2  +  2.9157xl08/V2 
The  military  thrust  will  be  computed  as  discussed  in  Chapter  2. 

T8,000ft  =  Tssl(^1/2) 

6  =  P8,000ft  /  Pssl 
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=  1571.88  2116.22 

=  0.742777 


G  = 


Ps.OOOft     Pssl 


-0.0018683  0.0023769 
=  0.786024 

T8,000ft  =  21,800  x  (0.742777/0.786024 12) 

=  18,2641b 

Fuel  flow; 

Fj  =  26,000-  0.64  x  Altitude 

=  26,000  -  0.64x8,000 

=  20,880 
F2  =  (Fj  +  333.3)  -  22.1  x  Drag  Index 

=  21,213.3  -22.1  x40 

=  20,329.3 
F3  =  F2  -  0.002  x  (Gross  weight  -  40,000) 

=  20,329.3  -  0.002  x  (50,000  -  40,000) 

=  20,309.3 

C  =  Military  thrust  /  Fuel  flow 
=  18,264/20,309.3 
=  0.8993 

Fuel  flow  =  Thrust  required/0.8993 

Table  17  and  Figure  3.6  show  the  results  of  calculation  and  deviation  from  hand 
book.  data. 

F.       VARIATION  NUMBER  5  -  WEIGHT,  DRAG  INDEX  AND  ALTITUDE 

We  have  seen  how  the  fuel  flow  was  changed  with  various  conditions.  Not  much 
time  is  needed  to  fmd  the  fuel  flow  for  conditions  which  can  be  read  directly  from  hand 
book,  but  most  of  the  conditions  do  not  correspond  to  the  published  conditions  which 
can  be  read  directly  from  the  hand  book  data  sheets.  Thus  to  find  the  fuel  flow  for 
specific  conditions  that  are  not  printed  in  the  hand  book,  one  must  spend  lots  of  time 
to  interpolate  each  item.  The  specific  fuel  flow  will  now  be  computed  by  both  methods. 
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Given  condition 

Gross  weight 

48,000  lbs 

Drag  Index 

27 

Altitude 

11,000  ft 

Velocity 

455  KTS 

The  procedure  to  compute  the  Kj  and  K2  is  the  same  as  in  the  previous  section. 
Consider  each  different  condition. 

1.  Gross  weight 

K7  =  K?  x  (48,000/40,000)2-7 

248,0001bs  240,0001bs 

=  2.0xl08x(1.2)2-7 

=  3.27xl08 

Ki  is  independent  of  changing  gross  weight. 

2.  Drag  index 

K{  =  0.0135  +  0.000132  x  Drag  Index 
=  0.0135  +  0.000132x27 
=  0.01706 

K2  is  independent  of  changing  Drag  Index. 

3.  Altitude 
1DI27 


Kl  "  Klni?7  X  (pl  1.000ft  /  Pssl)  ' 


=  0.01706  x  (0.0017007/0.0023769)1,2 


=  0.01142 

K2=  K248,0001b/(Pll'000ft/Pssl)1*2 
-  3.27xl08/(0.0017708/0.0023769)L2 

=  4.8893xl08 

Thrust  required  =  0.01142V2  +  4.8893xl08/V2 


=  0.01 142  x(455xl.69)2  +  4.8893xl08/(455xl.69)2 

=  7,572  lb 


Compute  the  factor  C: 
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Military  thrust  =  21,800  -  0.4  x  Altitude 
=  21,800-0.4x11,000 
=   17,400  lb 

Fuel  flow; 

Fj  =  20,840  -  0.37  x  (Altitude  -  8,000) 

=  20,840  -  0.37  x  (1 1,000  -  8,000) 

=  19,730 
F2  =  (Fj  +333.3)  -  22.1  x  Drag  Index 

=  20,063.3  -  22.1  x  27 

=  19,466.6 
F3  =  F2  -  0.002  x  (Gross  weight  -  40,000) 

=  19,466.6  -  0.002  x  (48,000  -  40,000) 

=  19,450.6 

C  =  Military  thrust  /  Fuel  flow 
=  17,400/19450.6 
=  0.89457 

Thus 

Fuel  flow  =  Thrust  required  /  C 
=  7,572/0.89457 
=  8,464  lbs/hr 

Now  the  fuel  flow  for  given  conditions  will  be  computed  by  interpolating  the 
hand  book  data.  Figure  3.7  shows  the  whole  diagram  to  compute  the  desired  fuel  flow 
of  given  conditions. 

Computations  for  the  steps  are: 

1)  Fl  =  (8873  -7988)x(27  -  20)/(40  -  20)  +  7988  =  8297.75 

2)  F2  =  (10324  -  9223)x(27  -  20)/(40  -  20)  +  9223  =  9608.35 

3)  F3  -  (455  -  440)x(9608.35  -  8297.75)/(480  -  440)  +  8297.75  =  8789.23 

4)  Gl  =  (7729  -  6962)x(27  -  20)/(40  -  20)  +  6962  =  7230.45 

5)  G2  =  (8840  -  7933)x(27  -20)/(40  -20)  +  7933  =  8250.45 

6)  G3  =  (455  -  440)x(8250.45  -  7230.45)/(480  -  440)  +  7230.75  =  7613.25 

7)  L4  =  (11000  -  8000)x(7613.25  -  8789.23)/(12000  -  8000)  +  8789.23  =  7907.24 

8)  HI  =  (9239  -  8366)x(27  -  20)/(40  -  20)  +  8366  =  8671.55 


80 


DIAGRAM  FOR  DESIRED  FUEL  FLOW 


GROSS  WEIGHT    40,000  LBS 

ALTITUDE   8,000  FT 

11,000  FT 

ALTITUDE   12,000  FT 

A/S 
KTAS 

DRAG  INDEX 

L4 

A/S 

KTAS 

DRAG  INDEX 

20 

27 

40 

20 

27 

40 

440 

7988 

Fl 

8873 

440 

6962 

Gl 

7729 

455 

F3 

455 

G3 

480 

9223 

F2 

10324 

480 

7933 

G2 

8840 

GROSS  WEIGHT    50,000  LBS 

ALTITUDE   8,000  FT 

11,000  FT 

ALTITUDE   12,000  FT 

A/S 
KTAS 

DRAG  INDEX 

M4 

A/S 
KTAS 

DRAG  INDEX 

20 

27 

40 

20 

27 

40 

440 

8366 

Hl 

9239 

440 

7447 

Jl 

8192 

455 

H3 

455 

^3 

480 

9522 

H2 

10688 

480     8315 

J2 

9216 

G/WT 

40,000  LBS 

48,000  LBS 

G/WT 

50,000  LBS 

A/S 
KTAS 

DRAG  INDEX 

N4 

A/S 
KTAS 

DRAG  INDEX 

27 

27 

455 

L4 

455 

M4 

Figure  3.7     Diagram  for  Desired  Fuel  Flow. 
81 


9)  H?  =  (10688  -  9522)x(27  -  20)/(40  -20)  +  9522  =  9930.1 

10)  H3  =  (455  -  440)x(9930.1  -  8671.55)/(480  -440)  +  8671.55  =  9143.51 

11)  Jl  =  (8192  -  7447)x(27  -  20)/(40  -  20)  +  7447  =  7707.75 

12)  J2  =  (9216  -  8315)x(27  -  20)/(40  -  20)  +  8315  =  8630.35 

13)  J3  =  (455  -  440)x(8630.35  -  7707.75)/(480  -  440)  +  7707.75  =  8053.73 

14)  M4    =    (11000   -    8000)x(8053.73    -    9143.51)/(12000   -    8000)    +    9143.51 
8326.174 

15)  N4    =    (48000   -   40000)x(8326.17   -    7907.24)/(50000-   40000)    +    7907.24 

Desired  fuel  flow  for  given  condition  is  8,242  lbs/hr. 
Thus  deviation  of  the  results 

(8,464  -  8,242)/8,242  x  100  =  2.7  % 
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IV.  CONCLUSION 

It  is  concluded  from  the  investigations  that  computing  the  fuel  flow  rate  from  the 
computer  program  results  in  less  deviations  and  is  simpler  than  interpolating  from 
NATOPS  FLIGHT  MANUAL  DATA.  The  total  lines  used  in  this  program  is  53  lines 
for  the  m:.in  frame  program  and  229  lines  for  the  hand-held  calculator  program,  and 
both  programs  are  inserted  in  Appendix. 

Though  the  deviation  of  fuel  flow  for  this  program  was  restricted  to  a  maximum 
of  10  percent,  for  most  conditions  deviations  are  distributed  below  5  percent.  To 
determine  the  fuel  flow  at  specific  conditions  which  can  not  be  read  directly  from  the 
flight  manual,  it  will  be  necessary  to  calculate  fifteen  steps  of  interpolations.  This  is 
very  time  consuming  and  the  probability  of  making  a  mistake  is  very  high  when 
calculating  so  many  interpolations. 

Thus  the  computer  program  to  compute  the  fuel  flow  rate  can  not  only  save  time 
for  many  conditions  which  must  be  changed  continuously  during  flight,  but  it  also 
reduces  the  probability  of  making  a  mistake  in  interpolations  steps. 

It  is  recommended  that  a  similar  computer  program  be  developed  for  specific 
flight  stages,  i.e.,  take-off,  climb,  or  landing,  using  the  same  concepts  that  were  used  in 
this  program. 
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APPENDIX  A 
COMPUTER  PROGRAM  (BASIC) 


00010  K2=2.0E+08  !  REM  xx  BASED  ON  S.S.L.  AND  GROSS  WEIGHT  40000LBS  XX 

00020  PRINT  'INPUT  ALTITUDE, GROSS  WEIGHT, DRAG  INDEX, AIR  SPEED, A/S  INC' 

00030  INPUT  A,W,I9,V1,VP 

00040  Kl=(0.0135+. 000132X19)  !  REM  XX  Kl  VARIES  WITH  DRAG  INDEX  xx 

00050  D4=1-6.875E-6*A  !  REM  XX  ICAO  SCALE  BASE  XX 

00060  T1=518.688XD4  !  REM  xx  ABSOLUTE  TEMP  (R)  BASED  ON  ALTITUDE  XX 

00070  D5=D4xx5.2561  !  REM  xx  PRESSURE  RATIO  BASED  ON  ALTITUDE  xx 

00080  A9=(l. 4*1714. 87*T1)  !  REM  xx  SPEED  OF  SOUND  (SQUARED)  xx 

00090  K9=(D5/D4)XX1 .2  !  REM  xx  ALTITUDE  CORRECTION  FOR  'Q»  xx 

00100  K1=K1XK9  !  REM  xx  PAXMXX2  ALTITUDE  CORRECTION  xx 

00110  K2=K2/K9  !  REM  xx  PAXMXX2  ALTITUDE  CORRECTION  xx 

00130  K2=K2X(W/40000)XX2.7  !REMxx  WEIGHT  OTHER  THAN  40000LB  CORRECTIONXX 

00160  PRINT  '  Kl  =  ' ■ K 1  'K2=,-K2 

00170  PRINT  TABU); ' VEL ( KTS ) ' ; TAB( 12) ; 'MACH' ; TAB( 20 ) ; 'PAR  DRAG'; 

00180  PRINT  TAB(34);'PRO  DRAG' ; TAB( 48 ); 'MACH  DRAG' ; TABC62 );' FUEL  FLOW' 

00190  V2  =  V1X1  .68894 

00200  M=V2/A9XX0.5 

00210  D8=l  !  REM  xx  CDO  MACH  CORRECTION  BASE  xx 

00220  IF   M>.8   THEN  480 

00230  D1=K1XV2XX2  !  REM  xx  PARASITE  DRAG  (SUBSONIC)  xx 

00240  D3=D1XD8  !  REM  xx  PARASITE  DRAG  (MACH  CDO)  xx 

00250  D2=K2/V2xx2  !  REM  xx  INDUCED  DRAG  xx 

00260  D0=D3+D2  !  REM  xx  MACH  DRAG  xx 

00270  D7=D1+D2  !  REM  XX  SUBSONIC  TOTAL  DRAG  xx 

00280  T=21800-AX0.4  !  REM  xx  MIL  THRUST  OTHER  THAN  S.S.L.  CORRECTION  xx 

00290  IF  A<=8000  THEN  320 

00300  F1=20840-(A-8000)X0.37!  REM  xx  F/FLOW  CORRECTTION  ABOVE  8000  FTxx 

00310  GO  TO  330 

00320  F1=26000-0.64000XA  !  REM  xx  F/FLOW  CORRECTION  BELOW  8000  FT  xx 

00330  IF  I9>20  THEN  360 

00340  F2=F1-16. 65X19  !  REM  xx  F/FLOW  CORRECTION  BELOW  DRAG  INDEX  20  xx 

00350  GO  TO  370 

00360  F2=(Fl+333.3)-22. 10X19  !  x  F/FLOW  CORRECTION  ABOVE  DRAG  INDEX  20  x 

00370  IF  W<=40000THEN  400 

00380  F3=F2-(W-40000)X0.002  !  REM  xx  F/FLOW  CORRECTION  WITH  WEIGHT  xx 

00390  GO  TO  410 

00400  F3=F2 

00410  C=T/F3  !  REM  xx  COEFFICIENT  FOR  FUEL  FLOW  xx 

00420  FF=D0/C  !  REM  xx  FUEL  FLOW  xx 

00430  PRINT  TAB(3);V1;TAB(8);M;TAB(18);D2;TAB(32);D3;TAB(46);D0; 

00440  PRINT  TAB(60);FF 

00450  V1=V1+VP 

00460  IF   Vl>610   THEN  510 

00470  GOTO  190 

0  048  0  D9=.55X(SIN((M-.8)/.8)XX3)/.0180 

00490  D8=l+D9 

00500  GOTO  230 

00510  PRINT  'C'^iC, 'T=»;T 

00520  STOP 

00530  END 
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APPENDIX  B 


COMPUTER  PROGRAM  (HAND-HELD  CALCULATOR) 


ei*LBL  "FF" 

51  PCL  63 

161  PCL  15       I 

CI 

+ 

281  RCL  23 

e:  i 

52  * 

162  RCL  13       ! 

c  2 

PCL  19 

262  -FF=- 

63  STO  24 

53  STO  69 

163  +          ] 

+ 

263  RRCL  X 

84  RCL  69 

54  RCL  64 

164  STO  17       1 

54 

STO  29 

264  RVIEH 

f=;  -flLT=> 

55  RCL  69 

165  PCL  96       1 

55< 

►LBL  "EE" 

265  STOP 

*i   PROMPT 

56  / 

166  .4          1 

56 

PCL  61 

266  "HEW  V?<1,9 

87  STO  68 

57  STO  64 

167  * 

57 

46869 

267  PROMPT 

63  RCL  61 

58  2.6  E63 

188  CHS 

58 

- 

268  X=8? 

69  "«=?" 

59  RCL  69 

169  21888        1 

[59 

.882 

269  STOP 

16  PROMPT 

66  / 

116  ♦ 

66 

* 

216  GTO  "VV 

11  STO  61 

61  STO  18 

111  STO  13       1 

61 

CHS 

21MLBL  -UK' 

12  RCL  62 

62  RCL  81 

112  8688        1 

62 

PCL  28 

212  RCL  12 

13  -DI=?' 

63  46868 

113  PCL  69 

63 

+ 

213  .8 

14  PROMPT 

64  / 

114  X<=Y?        1 

64 

STO  21 

214  - 

15  STO  62 

65  2.7 

115  GTO  -PP- 

65 

1/X 

215  .3 

16*LBL  "VV 

66  YtX 

116  rcl  ee-- 

66 

PCL  18 

216  / 

17  RCL  63 

67  RCL  19 

117  8699        ] 

67 

* 

217  EHTERt 

IS  -V<KTS>=?" 

68  * 

118  -          I 

63 

STO  22 

218  3 

19  PROMPT 

69  STO  18 

119  .37 

69 

1/X 

219  YtX 

26  STO  e3 

78  RCL  63 

126  *          1 

79  RCL  16 

226  SIN 

21  RCL  62 

71  1.68394 

121  CHS         I 

71 

* 

221  .55 

22  .666132 

72  * 

122  28848 

172 

STO  23 

222  » 

23  * 

73  STO  11 

123  + 

73 

PCL  63 

223  .618 

24  .6135 

74  RCL  67 

124  STO  19 

74 

FIX  8 

224  / 

25  ♦ 

75  S6RT 

125  GTO  -LL- 

1  J 

■V<KTS>= 

225  1 

26  STO  64 

76  !/X 

126»LBL  "PP" 

76 

RRCL  X 

226  + 

27-6.375  E-96 

77  RCL  11 

127  RCL  89 

77 

RVIEH 

227  STO  24 

23  RCL  69 

78  * 

12S  .64 

[78 

STOP 

223  GTO  'PR- 

29  * 

79  STO  12 

129  * 

79 

RCL  12 

229  .END. 

36  CHS 

86  .8 

136  CHS 

86 

FIX  3 

31  1 

81  X^=Y? 

13!  26966        1 

81 

•&RCH=- 

32  + 

82  GTO  "MH' 

132  +          1 

Oi. 

RRCL  X 

33  STO  65 

33»L6L  -RR- 

133  STO  19       1 

07 

RVIEH 

34  518.68 

84  RCL  11 

134*LBL  'LL*      1 

84 

STOP 

7C  * 

65  Xt2 

135  26          1 

85 

RCL  15 

36  STO  66 

86  RCL  64 

136  RCL  92       1 

86 

FIX  1 

37  1716.5 

87  * 

137  X>Y?        1 

87 

"D<I)=- 

■70  . 

•J  J  * 

88  STO  13 

138  GTO  -I'D*      1 

33 

RRCL  X 

39  1.4 

89  PCL  24 

139  16.65        1 

89  RVIEH 

46  * 

98  * 

149  * 

98 

STOP 

41  STO  67 

91  STO  14 

141  CHS          1 

91 

RCL  14 

42  PCL  65 

92  RCL  11 

142  PCL  19        I 

92 

•H<P)=- 

43  EHTERt 

93  XT2 

143  ♦           1 

93 

RRCL  X 

44  5.2561 

94  1/X 

144  STO  29        1 

94 

RVIEH 

45  YtX 

95  RCL  18 

145  GTO  -EE- 

95 

STOP 

46  STO  63 

96  * 

146»LBL  TiD'      1 

96 

RCL  16 

47  1116 

97  STO  15 

147  22.1 

97 

-D(H)=- 

48  Xt2 

98  PCL  14 

143  * 

98 

RRCL  X 

49  RCL  67 

99  ♦ 

149  CHS          I 

99 

RVIEH 

56  / 

166  STO  16 

156  333.3        t 

!88 

STOP 
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